1 



(N 
O 
O 
(N 

< 

(N 



Hadronic Spectral Functions above the QCD Phase Transition * 
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We extract the spectral functions in the scalar, pseudo-scalar, vector, and axial-vector channels above the 
deconfinement phase transition temperature (Tc) using the maximum entropy method (MEM). We use anisotropic 
lattices, 32^ x 32, 40, 54, 72, 80, and 96 (corresponding to T = 2.3Tc — > O.STc), with the renormalized anisotropy 
^ = 4.0 to have enough temporal data points to carry out the MEM analysis. Our result suggests that the spectral 
functions continue to possess non-trivial structures even above Tc and in addition that there is a qualitative change 
in the state of the deconfined matter between 1.5Tc and 2Tc. 
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1. INTRODUCTION 

The spectral fmictions (SPFs) of hadronic op- 
erators play an important role in QCD. The mod- 
ification of hadrons, which has been suggested, 
for example, by the dilepton production enhance- 
ment in high energy nuclear collisions observed at 
CERN SPS, can be formulated in terms of SPFs. 
Although there have been numerous theoretical 
attempts to understand the modification of SPFs 
at finite temperature and/or density Q, the ex- 
act nature of the hadronic modes in matter is not 
understood well. 

Recently, the maximum entropy method 
(MEM) has been used to extract SPFs from lat- 
tice QCD data for the first time by the present 
authors |^,^. The Euclidean 2-point functions in 
the temporal direction and the associated SPFs 
are related by the Laplace transform. On the 
lattice, only a finite number of data points with 
statistical noise are available in the temporal di- 
rection. Therefore, direct Laplace inversion from 
lattice data to SPF is an ill-posed problem. MEM 
is a method to evade such difficulty on the basis 
of the Bayes' theorem in the theory of statistical 
inference Hi. 



2. HOW MANY DATA POINTS 
NECESSARY FOR MEM ? 
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As we have discussed in detail in 0, the re- 
sult of MEM depends strongly on the number of 
temporal data points. The more data points are 
used, the closer the MEM result is to the true 
SPF. Moreover, there is the minimum number of 
data points TV to perform the reliable MEM anal- 
ysis. 

In order to find this minimum at T = on 
the lattice, we have carried out the following anal- 
ysis [Q. First we calculated hadronic correlators 
with the quenched approximation on an isotropic 
40^ X 30 lattice at ^ = 6.47 The number 

of gauge configurations is 160. Then, MEM anal- 
ysis has been done for the vector channel by us- 
ing N data points (at r = 2, • ■ • , N/2 + 1 and 
31 — N/2, • ■ • , 30) out of the total 30 data points. 
The result is shown in Fig. |^. The figure clearly 
shows that, SPF changes considerably as N and 
that at least about 30 data points are necessary 
for the convergence of the result. This minimum 
N would depend on f3, and on whether one em- 
ploys anisotropic lattices, improved actions, and 
so on. Nevertheless, we use this number as a prac- 
tical guide in the following. 
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Figure 1. MEM results in the vector channel with 
10 different iV's. 
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Figure 2. Polyakov loop expectation value as a 
function of Nt . 



3. ANISOTROPIC LATTICE AT T ^ 

At finite temperature (T) , because of the rela- 
tion T = l/Nrttr {Nr and a,- being the temporal 
lattice size and spacing, respectively), less and 
less data points are available as T increases if Qr 
is fixed. Requiring that N-r > 30 holds even at 
the highest T, e.g. 2.5Tc, we are inevitably led to 
use an anisotropic lattice. 

We have used the bare anisotropy = 3.5 and 
/3 = 7.0 with the naive plaquette action. For 
the quark part, the standard Wilson action with 
the quenched approximation is used. The cor- 
responding renormalized anisotropy is ^ = 4.0 
1^ and the lattice spacing is Ct = 1/4 • Oo- = 
9.75 X 10"'^ fm. Simulations are done on 32^ x Nr 
lattices with = 32, 40, 54, 72, 80, and 96. This 
corresponds to T ~ 2.3re, l.QT^, lAT^, IMTc, 
0.93Tc, and 0.78rc, respectively. More than 100 
configurations are generated for each N^. The 
details of the lattice parameters used in the cal- 
culation will be given in j|] . 

We have studied 2-point correlation functions 
in the scalar (S), pseudo-scalar (PS), vector (V), 
and axial-vector (AV) channels. The default 
models, motivated by perturbative QCD, in each 
channel are m = 0.60, 1.15, 0.40, and 0.35, 
respectively. The lattice data points used in 
our MEM analysis are r — 5, • • • , 21 and Nr — 
19, • • ■ , iV^ - 3 for Nr = 54, 72, 80, and 96. For 



Nr = 32 and 40, r = 5,---,7Vt- - 3 are used. 
When the source and sink are closer than ^Ur in 
the temporal direction, artifact could arise due 
to unphysical modes at w > n/S^ar- Therefore, 
we leave out the data points near the edges, i.e., 
T = I,- - ■ ,4 and Nr - 2, - ■■ ,Nr. 

The number of lattice data points used in the 
MEM analysis is fixed almost the same at each 
Nr except at Nr = 32 to keep the resolution in 
the MEM analysis unchanged. The rest of the 
MEM procedures such as the averaging over a 
are basically the same as described in and will 
be presented in detail in 

Fig. H shows the Polyakov loop expectation 
value |(L)| at each Nr- It starts to deviate from 
zero around Nr = 80(72), which corresponds to 
Tc — 253(281) MeV, and is comparable to the real 
value 271 ± 2 MeV. We have checked that the 
Polyakov loop susceptibility also shows a peak in 
this region. In the following, we use the number 
Tc = 271 MeV. 

4. RESULTS OF MEM 

In Fig. |, we show the results of the MEM 
analysis for SPF p(w) in the S, PS, V, and AV 
channels on 32^ x 54 lattice (T ~ 1.4Tc). In this 
calculation, we have used the hopping parameter 
corresponding to m^/m^ ~ 0.7 at T = 0. 

If the quark-gluon plasma is such an object as 
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described as a free gas of massive quarks, anti- 
quarks and gluons, p{u)) would show a flat be- 
havior with a smooth rise from zero above the 
qq threshold as a function of uj. To the contrary, 
Fig. I shows that SPFs possess nontrivial struc- 
tures even above Tc in all four channels. There 
are a sharp peak at about u) = 2 GeV and two 
bumps above the peak. Also SPFs in all channels 
almost degenerate, which shows a strong signa- 
ture of the restoration of chiral symmetry. The 
peak and bump structures, and the suppression 
of SPFs below uj — 2 GeV are statistically sig- 
nificant according to our error analysis. We have 
used Nr twice as large as that employed in 0. 
As shown in Sec. 2, this is necessary for reliable 
MEM analysis. 

In Fig. 1^, we show the results at Nr = 40, i.e., 
T ~ 1.9Tc. There is an apparent peak around 
w = 0. A possible explanation of the peak is the 
effect of Landau damping ||]. However, this peak 
is, at the moment, not statistically significant yet. 
The peaks around oj ~ 4.5 GeV are significantly 
broader than those around 2 GeV at T ~ l.OTc- 

We have carried out an error analysis for the 
average of the SPFs in small regions around the 
peaks and found that the broadening is indeed 
statistically significant. SPFs in the case of Nr — 
32 (T ~ 2.3rc) have similar structures as those at 
T ~ 1.9Tc except the shift of the peak positions. 
This suggests a possibility of a qualitative change 
in the state of the deconfined phase between l.bTc 



and 2Tc. 
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